AO  AOS  1 1 08 


'#  ft. 

DNA  3929F 

DNA  PROJECT  609  RADAR: 

AURORAL  BACKSCATTER  MEASUREMENTS 

Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park,  California  94025 


February  1976 


Final  Report  for  Period  31  January  1974—30  September  1975 


CONTRACT  No.  DNA  001-72-C-0179 


Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  D.  C.  20305 


UNCLASSIFIED 

PHTv  4 ' a • ..  F Th  PACI  Whmn  D«<  I 

/J I"’  REPORT  DOCUMENTATION  PAGE 

1 Ot  u T,  /-./-ve/t  Ar  r r«.cir»u  3 Rf  nT’U  *T  AlOC  NUMBfS 

^ DNa)  p929F*~7 / 1 I ^7  y 


R t n ^ H 


I 4 T I TlF  unit  '..Mill. 


4TS  PERIOD  COVERED 


/ DNA. PROJECT  60 

^ ^ASUREMENTS  , 


609  RADAR:  AURORAL  BACKSCATTER  I Final  Rep«*t,f«»  Taui 

/ I 31  Jan  74  — 3/8  Sep  75 

SRI  Project  2001  V— 


I 7 A T H R 


I " R.  T. /Tsunoda . 
( R.  I . Presnell 

T.  N.  C./toanq 


>N  N AM  t AND  ADORE  SS 


Stanrord  Research  Institute  I ( *i 

333  Ravenswood  Avenue  ' v NWET-SW***#!* 

Menlo  Park,  California  94025  |\  L25AAXHj— — 

1’  CONTROLLING  OFFICE  NAME  AND  ADDRESS  .x—  I 

Director  j | ‘ ( 

Defense  Nuclear  Agency  “ •< 

Washington,  D.C.  20  305  54  ^ 

~il  m.gniT  1R|«|J  AVER  uy  NAME  0 4DORESS.il  Ittlrrvnl  li.m.  I .nlr.i  Hint  nil"  '■  IT?-  ,1  uR|Tv 

_ 

I y t , “ UNCLASSIFIED 


0 CONTRA  T OR  GRANT  4 jUBI  4 


L DNA  001-72-C-0179 


10  PROGRAM  E.EMEnT  PROJECT  Tag* 
- 4AEAA  4 |B»:  UNIT  N IMRERS 


NWET-fc 

L25AAX 


1?  / 1 


IS*  DECLASSIFICATION  DOWNGRADING 


16  Distribution  STATEMENT  (of  rh»s  Report  > 


Approved  for  public  release;  distribution  unlimited. 


’•  ' T OieuTJON  STATEMENT  'ol  thr  Abstract  entered  in  Block  20.  tl  different  from  Report) 


o * s 

r4' % 

€.:^o 


10  SUPPL EMEN T ARY  NOTES 


This  work  sponsored  by  the  Defense  Nuclear  Agency  under  Subtask 
L25AAXHX6 31-07. 


!*»  kf  v WORDS  i < nntinue  .n  reverse  side  if  necessary  and  identify  by  block  number  i 


Ionospheric  Physics 
Auroral  Radar  Clutter 
Phased-array  Radar 
Auroral  Ionosphere 


Ionospheric  Electric  Fields 
Ionospheric  Motions 
Electrojet  Currents 
Auroral  Storms 


i . i. ; ■ . rni.-  ,.n  rr .....  4,.  ft.rf  Ki*nn  h At  blrult  m/nthrr) 

In  this  report summarizejthe  research  results  obtained  from  auroral 
backscattT?f  studies  using  the  DNA  609  radar  system  located  at  Homer,  Alaska. 
The  primary  purpose  of  the  b^ckscatter  studies  was  to  improve -«*•»“  under- 
standing of  the  auroral-clutter-producing  mechanism.  Most  of  the  results 
are  based  on  coordinated  experiments  utilizing  the  398-MIIz  phased-array 
radar  at  Homer,  the  SRI/DNA  incoherent-scatter  radar  at  Chatanika,  and 
other  complementary  auroral  measurement  techniques.  We  found  that  the 
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20.  ABSTRACT  (Continued) 

problem  can  be  most  effectively  approached  by  seeking  both  microscale  and 
macroscale  relationships--i . e. , the  relationships  of  local  plasma  para- 
meters within  the  radar  scattering  volume  producing  auroral  clutter,  and 
the  relationships  of  the  large-scale  auroral  and  magnetospheric  conditions 
to  the  spatial  distribution  of  auroral  clutter,  respectively.  The  results 
include  the  following: 


(1)  Diffuse  auroral  clutter  is  apparently  not  produced  by  primary 
plasma  waves  generated  by  the  Buneman-Farley  two-stream 
instability. 


(2)  The  mechanism  does  not,  however,  appear  to  be  unlike  the  quasi- 
linear  gradient-drift  instability  model  proposed  by  Sudan  et  al . 
(1973),  which  can  be  driven  by  the  primary  two-stream  waves. 


(3)  There  is  a threshold  electric  field  strength  of  30  mV/m  associated 
with  the  occurrence  of  398-MHz  diffuse  auroral  clutter. 


(4)  The  diffuse  auroral  clutter  is  closely  associated  with  the 
auroral  electrojets. 

(5)  The  evening  diffuse  auroral  clutter  is  closely  associated  with 
downward  field-aligned  currents. 

(6)  Discrete  auroral  clutter  is  closely  associated  with  the 
discrete  visual  aurora. 

In  addition  to  the  above  rfehlts,  we  present  (1)  a brief  data  summary 
taken  in  support  of  the  DNA  ICECA?  '74A  rocket  program;  (2)  a description 
of  the  near-on-line  Doppler  capability  that  was  added  to  the  phased-array- 
radar  system;  (3)  a discussion  of  A\field  trip  to  search  for  plasma-line 
backscatter  associated  with  auroral  asftivity;  and  (4)  the  current  status 
of  a theoretical  study  considering  the '^on-cyclotron  instability  as  a 
contributor  to  the  production  of  aurorai-Nclutter . Finally,  we  discuss  the 
need  for  future  studies  of  the  characteristics  of  F- region  auroral  clutter, 
a phenomenon  that  has  received  little  or  ncS^ittention  in  the  past.  This 
type  of  auroral  clutter  has  recently  received  theoretical  support  in  the 
form  of  a plasma  microinstability  model. 
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I INTRODUCTION 


DNA  Project  609  originated  in  the  fall  of  1961.  At  that  time,  the 
design  and  construction  of  clutter  radars  was  started  in  preparation  for 
the  FISH  BOWL  nuclear  test  series  the  following  year.  Following  that 
test  series,  the  radars  were  moved  to  Homer,  Alaska  to  conduct  auroral 
clutter  studies.  The  radar  system  at  that  time  consisted  of  six  oper- 
ational 1 requencies — 50,  139,  398,  850,  1210,  and  3000  MHz — all  of  which 
were  transmitted  and  received  through  a common  60-ft  parabolic  dish 

antenna.  The  research  results  of  past  studies  have  been  published  in  a 

1,2* 

number  of  Project  609  final  reports. 

In  1970,  the  design  of  the  398-MHz  phased-array  antenna  system  was 
started.  The  array  was  first  assembled  and  tested  at  Stanford,  California 
during  1972.  The  array  was  then  disassembled,  shipped,  reassembled,  and 
integrated  with  the  398-MHz  radar  system  at  Homer,  Alaska.  The  first 
operation  of  the  phased-array  radar  in  an  auroral  environment  was 
conducted  during  the  spring  of  1973. 

It  became  immediately  obvious  from  the  first  operation  that  the 

phased-array  radar  provided  a quantum  jump  forward  in  our  ability  to 

unravel  the  lingering  problems  concerning  the  interpretation  of  aurora  , 

► 

clutter-producing  mechanisms.  By  the  radar's  rapid-scan  capability  and 
on-line  computer  support,  the  spatial  distribution  of  auroral  clutter, 
with  minimal  temporal/spatial  ambiguities,  could  be  mapped  over  a large 

* 

References  are  listed  at  the  end  of  this  report. 
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geographical  area.  From  these  maps,  true  clutter  boundaries  (and  their 
dynamics)  could  be  discerned  that  were  related  to  auroral  and  magneto- 
spheric  processes,  and  that  were  not  due  to  radar  characteristics. 

Equally  important  was  the  installation  of  the  Chatanika  incoherent- 
3 

scatter  radar  beneath  the  central  region  scanned  by  the  Homer  radar. 

With  the  Chatanika  radar,  plasma  parameters  such  as  the  electric  field, 
electron  density,  and  current  density  could  be  measured  within  the  Horner 
radar  scattering  volume.  With  this  combination,  complemented  by  other 
auroral  measurement  techniques  (e.g.,  all-sky  cameras,  magnetometers,  and 
rocket-borne  and  satellite-borne  instrumentation),  great  strides  have 
been  made  toward  an  understanding  of  the  auroral-clutter-producing 
mechanism . 

The  first  research  results  using  the  1973  data  were  presented  in 

1 

the  most  recent  Project  609  final  report.  Since  then,  several  more 
field  experiments  have  been  conducted.  Data  analysis  was  continued  on 
both  the  1973  data  and  those  collected  during  the  subsequent  field 
experiments.  The  research  results  covering  the  time  period  from 
31  January  1971  to  30  September  1975  are  the  principal  subject  of  this 
report . Details  of  the  various  studies  conducted  under  this  project 
have  been  published  in  the  following  journal  articles: 

• Tsunoda,  K.  T. , R.  I.  Presnell,  and  R.  L.  beadabrand,  "Radar 
Auroral-Echo  Characteristics  as  Seen  by  a 3 98 -MHz  Phased- 
Array  Radar  Operated  at  Homer,  Alaska,”  ■)  . Geophys.  Res., 

79,  p.  1709,  1971. 

• Wang,  T.  N.  C.  and  R.  T.  Tsunoda,  "On  a Crossed  Field  Two- 
.Stream  Plasma  Instability  in  the  Auroral  Plasma,"  J . Geophys. 
Res. , 80,  p.  2172,  1975. 

• Hunsucker,  R.  D. , G.  J.  Romick,  W.  L.  Ecklund,  R.  A.  Greenwald, 
B.  B.  Balsley,  and  R.  T,  Tsunoda,  "Structure  and  Dynamics  of 
Ionization  and  Auroral  luminosity  During  the  Auroral  Events 

of  March  16,  1972,  near  Chatanika,  Alaska,"  Radio  Sci . , 10, 
p.  813,  1975. 
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• Tsunoda,  f(.  T.,  "Electric  Field  Measurements  Above  a Radar 
Scattering  Volume  Producing  'Diffuse*  Auroral  Echoes, " 

J.  Geophys.  lies.,  80,  p.  4297,  1975. 

• Tsunoda,  Ft.  T.,  "Doppler  Velocity  Maps  of  the  Difluse  Kadar 
Aurora,"  .J  . Geophys.  lies.,  81,  in  press,  1976. 

• Tsunoda,  It.  T.  and  ft.  I.  Presnell,  On  a Ttireshold  Electric 
Field  Associated  with  the  398-MHz  Diffuse  Kadar  Aurora," 

J,  Geo  pity  s,  lies.,  81 , p.  88,  1976. 

• Tsunoda,  It.  T.,  R.  1.  Presnell,  and  T.  A.  Potemra,  "The  Spatial 
Relationship  Between  the  Evening  Radar  Aurora  and  Field- 
Aligned  Currents,"  J . Geophys . Res . , 81 , in  press,  1976. 

• Tsunoda,  It.  T.,  R.  I.  Presnell,  Y.  Kamide,  and  S.  -I.  Akasofu, 
"Relationship  of  Radar  Aurora,  Visual  Aurora,  and  Auroral 
Electrojets  in  the  Evening  Sector,”  submitted  to  J,  Geophys.  Res., 
1976. 

Papers  in  preparation  include: 

• Wang,  T.  N.  C.  and  R.  T.  Tsunoda,  "The  Effects  of  Field- 
Aligned  Gradients  on  the  Two-Stream  Instability  in  an 
Auroral  Plasma." 

• Brekke,  A.,  R.  T.  Tsunoda,  and  M.  J.  Baron,  "On  the  Azimuthal 
Dependence  of  Slant  Range  Doppler  Shift.” 

Because  most  of  the  research  results  have  appeared  or  will  appear  in  the 
open  literature,  we  present  only  a summary  of  the  results.  The  results 
of  the  data  analysis  are  presented  in  Section  II.  A description  of 
field  experiments  (in  particular,  support  of  the  I)NA  ICECAP  '74  rocket 
program)  and  any  additions,  improvements,  or  modifications  to  the  radar 
system  are  described  in  Section  III.  The  results  of  theoretical  studies 
evaluating  the  ion-cyclotron  instability  as  a possible  source  of  auroral 
clutter  are  presented  in  Section  IV.  A discussion  of  the  results  is 
presented  in  Section  V.  And  finally,  conclusions  and  recommendations 
for  further  studies  are  presented  in  Section  VI. 
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II  DATA  ANALYSIS 


In  the  study  of  auroral-clutter-producing  mechanisms,  we  have  found 
that  the  problem  is  most  effectively  approached  by  seeking  both  "micro- 
scale" and  "macroscale"  relationships.  By  microscale,  we  refer  to  the 
local  conditions  within  the  radar  scattering  volume  that  seem  to  deter- 
mine the  behavior  of  auroral  clutter.  In  particular,  we  were  most  inter- 
ested in  plasma  parameters  such  as  the  electric  field,  current  density, 
electron  density,  etc.  within  the  radar  scattering  volume,  and  their 
relationship  to  auroral  clutter  characteristics.  And  by  macroscale, 
we  refer  to  the  large-scale  auroral  and  magnetospheric  conditions  that 
appear  to  be  associated  with  the  spatial  distribution  of  auroral  clutter. 
In  particular,  we  were  most  interested  in  relating  the  spatial  distri- 
bution of  the  convection  electric  field,  particle  precipitation,  auroral 
electro.jets , field-aligned  currents,  visual  aurora,  etc.,  to  the  spatial 
distribution  of  auroral  clutter.  Furthermore,  by  comparing  the  macro- 
scale relationships  with  microscale  relationships,  we  were  able  to 
derive  a consistent  phenomenology  from  which  we  can  more  intelligently 
interpret  the  auroral  clutter  characteristics,  particularly  under  complex 
and  dynamic  conditions. 

A.  Microscale  Relationships 

Comprehensive  microscale  studies  of  auroral  clutter  characteristics 

were  made  possible  by  the  installation  of  the  incoherent-scatter  radar 

3 

at  Chatanika,  Alaska.  The  locations  of  the  Homer  and  Chatanika  radars 
are  shown  on  a map  of  Alaska  in  Figure  1.  Chatanika  is  situated  directly 
beneath  the  central  region  scanned  by  the  Homer  radar  (shown  by  a pie- 
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REGION  SCANNED  BY 


kZIMUTM  - DISTANCE  PROJECTION 
FROM  HOMER 


shaped  sector).  Consequently,  the  incoherent-scatter  radar  can  be  used 
to  monitor  the  plasma  parameters  within  and/or  in  the  vicinity  of  the 
Homer  radar  scattering  volume. 


The  first  coordinated  study  using  the  two  radars  was  conducted  in 

5 

1972  and  the  results  were  published  in  a Chatanika  radar  final  report 

6,7 

and  subsequently  as  two  .journal  papers.  The  approach  taken  in  that 

study  was  to  compare  the  auroral  clutter  characteristics  with  the  electric 
field  vector  and  neutral  wind  determined  from  Chatanika  radar  measure- 

8,9  „ 

merits.  The  lirst  significant  result  was  the  finding  that  diffuse 

4,10 

auroral  clutter  was  not  produced  by  primary  plasma  waves  generated 

11,12 

by  the  Buneman-Farley  two-stream  instability.  However,  the  electric 

field  strength  was  found  to  be  large  enough  so  that  the  two-stream 

instability  was  operative.  In  other  words,  while  primary  plasma  waves 

were  present,  they  were  not  responsible  for  the  backscatter  observed 

over  the  azimuth  sector  scanned  by  the  Homer  radar,  From  the  analysis 

of  the  Doppler  spectra  of  auroral  clutter,  we  concluded  that  diffuse 

auroral  clutter  may  be  associated  with  secondary  plasma  waves  generated 

by  the  same  instability.  An  alternative  interpretation  is  that  these 

13 

secondary  waves  could  be  produced  by  the  gradient-drift  instability 

driven  by  the  primary  two-stream  plasma  waves.  The  interaction  between 

the  two  instabilities,  which  is  accounted  for  in  the  quasi-linear  plasma 

13 

instability  model,  was  first  proposed  by  Sudan  et  al.  and  applied  to 

14 

the  auroral  geometry  by  Greenwald.  The  interaction  was  suggested  by 

the  observed  Doppler  velocity  and  its  apparent  dependence  on  the  angle 

made  between  the  radar  beam  and  the  direction  of  current  flow.  Finally, 

6 

the  analysis  showed  that  the  occurrence  of  auroral  clutter  at  398  MHz 
might  be  associated  with  a threshold  electric  field  strength  of  approxi- 
mately 30  mV/m. 

The  possible  relationship  between  ionospheric  electric  field  strength 
and  the  occurrence  of  diffuse  auroral  clutter  was  confirmed  and  presented 
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in  the  previous  Project  609  Final  Report  and  in  a paper  by  Tsunoda  and 
15 

Presnell.  With  the  398-MHz  phased-array  radar,  dill  use  auroral  clutter 
was  lound  to  occur  over  Chatanika  whenever  the  measured  ionospheric 
electric  field  strength  exceeded  a nominal  value  of  30  mV/m.  When  the 
electric  field  strength  was  less  than  that  value,  auroral  clutter  was 
not  observed,  regardless  of  the  electron  density.  The  threshold  value 
is  consistent  with  the  activation  of  the  two-stream  instability  as  the 
primary  driving  mechanism. 

B.  Macroscale  Utlal lonship- 

The  first  results  of  this  kind  were  obtained  through  analysis  of 

4 

the  Spring  1973  data  and  reported  in  the  previous  linal  report  and  in  a 

10 

paper  on  the  first  phased-array  radar  results.  IVe  showed  that  auroral 

clutter  displayed  a consistent  morphology  based  on  spatial  distribution 

of  the  ionospheric  electric  field  and  the  E-region  electron  density.  The 

requirement  of  both  an  electric  field  and  an  electron  density  clarified 

the  hitherto  puzzling  relationship  between  auroral  clutter  and  the  visual 
1 

aurora.  The  diffuse  auroral  clutter  that  occasionally  extends  over  sever- 
al hundreds  of  kilometers  in  latitude  was  found  in  the  evening  sector  to 
be  located  equatorward  of  the  most  equatorward  visual  arc.  This  result 
is  apparently  due  to  (1)  subthreshold  electric  field  strengths  associated 
with  regions  of  intense  precipitation,  and/or  (2)  the  absence  of  signi- 
ficant particle  precipitation  in  regions  poleward  of  the  diffuse  auroral 
clutter  and  distant  from  the  locations  of  the  visual  arcs. 

16 

The  second  interpretation  is  supported  in  another  study  comparing 
the  spatial  distributions  of  auroral  clutter,  auroral  electrojets  and 
visual  aurora.  The  eastward  electrojet  (as  determined  by  a meridian  chain 
of  magnetometers)  was  found  to  be  collocated  with  the  diffuse  auroral 
clutter.  In  fact,  the  boundaries  of  diffuse  auroral  clutter  could  be 
used  to  define  boundaries  of  the  eastward  electrojet.  These  results  are 


9 


17 , 18 

consistent  with  the  findings  of  Greenwald  et  al . for  50  MHz,  that  the 

range-integrated  backscatter  amplitude  is  often  linearly  proportional  to 


the  horizontal  perturbation  component  \ (All)  + (AD)*'  of  the  geomagnetic 
field.  The  occurrence  of  the  eastward  electrojet  (and  diffuse  auroral 
clutter)  in  a region  equatorward  of  the  discrete  visual  aurora  is  con- 
sistent with  the  presence  of  large  electric  fields  in  regions  of  diffuse 

19 

particle  precipitation. 

In  the  morning  sector,  both  diffuse  and  structured  complex  auroral 

clutter  were  associated  with  the  regions  containing  visual  aurora.  The 

difference  in  the  spatial  correlation  of  auroral  clutter  with  visual 

aurora  implies  that  the  electric  field  strength  is  not  depressed  or  small 

within  visual  forms  that  occur  in  the  morning  sector.  Auroral  arc  models 

that  allow  polarization  fields  to  exist  within  arcs  and  in  which  the 

westward  electric  field  drives  a Pedersen  current  appear  to  be  consistent 
20 

with  the  data. 
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In  another  study,  we  found  that  the  diffuse  auroral  clutter  in 
the  evening  sector,  which  was  found  to  be  collocated  with  the  eastward 
electrojet,  was  also  collocated  with  the  region  containing  downward  field- 
aligned  currents.  The  field-aligned  currents  were  detected  by  a vector 
magnetometer  onboard  the  Navy/APL  TRIAD  satellite.  The  collocation  of 
field-aligned  currents  with  auroral  clutter  suggests  that  plasma  insta- 
bilities driven  by  such  currents  must  be  considered  as  possible  contri- 
butors to  the  production  of  auroral  clutter.  For  this  reason,  a 
theoretical  study  of  the  ion-cyclotron  instability  was  initiated.  The 
first  results  of  that  study  are  presented  in  Section  IV. 

The  relationships  discussed  above  can  be  used  to  interpret  what  is 

now  believed  to  be  a typical  evening-sector  case  under  moderately 

o 

disturbed  conditions.  An  example  is  shown  in  Figure  2.  The  40  azimuth 
sector  typically  scanned  by  the  Homer  radar  is  seen  to  contain  the  two 
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basic  types  of  auroral  clutter,  the  diffuse  band  (DB),  and  the  discrete 

arc  (DA).*’*0  Visual  arcs  scaled  from  a 11 -sky -camera  (ASC)  photographs 

taken  at  Fort  Yukon  and  Chatanika  are  shown  as  solid  lines.  (A  third, 

faint  visual  arc  seen  in  the  upper  half  of  the  Fort  Yukon  ASC  photograph 

was  not  scaled  because  there  was  no  corresponding  auroral  clutter  feature.) 

In  addition,  the  poleward  boundary  of  what  is  believed  to  be  the  diffuse 
22,23 

aurora  has  also  been  scaled  and  shown  as  a dashed  line. 

The  I)B  is  seen  to  extend  from  just  north  of  Talkeetna  (and  limited 
by  magnetic  aspect  sensitivity)  to  almost  directly  over  Fort  Yukon.  The 
DB  poleward  boundary  is  assumed,  in  this  case,  to  correspond  to  the 
poleward  boundary  of  the  weak  diffuse  echoes  that  is  coincident  with  the 
poleward  edge  of  the  diffuse  visual  aurora.  If  this  is  the  case,  the 
two  visual  arcs  and  the  DA,  which  is  situated  between  them,  must  be 
imbedded  within  the  DB. 

The  morphology  is  consistent  with  that  found  from  scanning  photo- 

22,23 

meter  data  taken  from  the  ISIS-2  satellite.  The  discrete  auroral 

arcs  are  imbedded  within  but  near  the  poleward  boundary  of  the  region  of 
diffuse  particle  precipitation.  The  I)B  corresponds  to  a region  of 
enhanced  electric  field  strength,  diffuse  particle  precipitation,  down- 
ward field-aligned  current  flow,  and  the  eastward  electrojet.  The  I)B 
is  modulated  near  its  poleward  boundary  due  to  a corresponding  modulation 
in  the  electric  field  produced  by  intense  localized  particle  precipitation 
associated  with  the  visual  arcs.  The  DB  poleward  boundary  is  therefore 
produced  by  the  absence  of  significant  particle  precipitation  poleward  of 
the  diffuse  aurora.  Upward  field-aligned  currents  were  found  to  be 
associated  with  the  discrete  visual  arcs.16’2** 
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Ill  FIEIJ)  EXPERIMENTS 

A.  ICECAP  ’71 

Two  field  trips  were  made,  to  participate  in  the  DNA  ICECAP  '74A 
and  '71U  rocket  programs.  The  Homer  radar  was  operated  between  Feb- 
ruary 13  and  March  2 during  the  1 7 1A  rocket  launches  from  the  Poker  Flat 
rocket  range  and  between  April  5 and  19,  the  time  period  of  the  ’7113 
rocket  program.  However,  due  to  on-line  computer  difficulties,  no  auroral 
clutter  data  were  collected  during  the  times  of  the  ICECAP  '7113  rocket 
launches.  Approximately  46  hours  of  clutter-map  data  were  collected 
during  the  two  field  operations.  In  the  following  subsections,  we 
present  a brief  summary  of  the  Homer  radar  data  collected  during  the  time 
periods  around  the  ICECAP  '74A  rocket  launches. 

1 . 14  February  1971 — Black  Brant  and  U te  Tomahawk 

A Black  Brant  rocket  (A18. 006-4)  carrying  instrumentation  to 

measure  IK  spectra  during  quiet  conditions  was  launched  from  Poker  Flat 

at  0607  UT . Thirty-seven  minutes  later,  a Ute-Tomahawk  rocket  was  also 

launched,  carrying  instrumentation  to  determine  the  neutral  atmospheric 

density  and  temperature,  also  during  quiet  conditions.  The  geomagnetic 

25 

and  auroral  conditions  during  this  time  period  were  relatively  quiet, 
with  the  College  H-component  only  about  10V  above  the  quiet-time  level. 

Any  visual  auroral  activity  that  may  have  been  over  Chatanika  was  obscured 
by  clouds. 

• Auroral  clutter  was  lirst  seen  around  0555  UT  approximately 

300  km  northeast  of  Chatanika,  The  echoes  were  relatively  limited  in 
both  latitudinal  and  longitudinal  extent,  in  contrast  to  the  broad 
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diffuse-band  (DU)  signatures  that  are  often  seen  during  this  time 

period  under  more  active  auroral  conditions  (see  Figure  2).  By  the 
time  of  the  rocket  launch,  the  1)B  had  moved  within  200  km  north  ol 
Chatanika  (anil  directly  over  Fort  Yukon),  but  had  also  weakened.  An 
auroral  clutter  map  taken  around  two  minutes  after  launch  is  shown  in 
Figure  3.  The  DB  is  seen  to  be  patchy  but  can  be  shown  to  be  more  or 
less  aligned  along  a constant  magnetic  1,-shell.  The  DB  can  be  seen,  by 
exercising  a little  imagination,  to  be  bifurcated.  That  is,  there 
appears  to  be  a region  within  the  DB  where  no  auroral  clutter  is  observed. 
That  region  also  appears  to  be  aligned  along  an  1,-shell.  The  auroral 
clutter  continued  to  weaken  and  was  essentially  gone  by  0606  UT. 


Peak  electron  densities  measured  by  the  Chatanika  radar  were 
•1  3 

typically  less  than  5 > 10  el/cm  in  the  E region  and  no  higher  than 
5 3 25 

10  el/cm  in  the  F layer.  Since  the  returned  incoherent-scatter  power 
was  so  small,  reliable  estimates  of  the  electric  field  could  not  be  made. 

15 

Tsunoda  and  Presnell  showed  that  in  addition  to  a threshold 

electric  field  strength  of  30  mV/m,  there  was  a minimum  electron  density 

below  which  diffuse  auroral  clutter  also  could  not  be  detected.  That 

4 , 3 

value  was  a few  parts  times  10  el/cm  . While  it  was  not  established 
whether  the  minimum  electron  density  represented  a true  "threshold"  in 
the  electron  density  or  whether  it  corresponded  with  the  sensitivity 
threshold  of  the  radar,  the  low  E-region  electron  density  measured  by 
the  Chatanika  radar  implies  that  the  electric  field  strength  could  have 
been  greater  than  30  mV/m.  If  this  is  the  case,  the  presence  of  bifur- 
cated DB  over  Fort  Yukon  may  be  due  to  the  presence  of  an  auroral  arc  in 
the  region  of  no  radar  echo.  The  auroral  arc  corresponds  in  location  to 
a region  of  localized  particle  precipitation,  and  hence,  enhanced  electron 
density.  The  absence  of  auroral  clutter  within  the  bifurcated  region 

could  lie  explained  by  the  depression  of  the  electric  field  (in  the  evening 

15 

sector)  by  the  intense  ionization.  This  question  can  be  answered  by 
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referring  t « > optical  (cither  all-sky  or  photomt.ar)  data  taken  from  Fort 

16 

\ukon.  Based  on  the  findings  presented  in  Tsunoda  et  al.,  that  visual 

arc,  if  present,  is  probably  located  near  the  lower  latitudinal  boundary 

26 

of  the  Harang  discontinuity. 

The  second  rocket  was  launched  at  064-1  IT  during  quiet  magnetic 
conditions  that  persisted  from  the  first  rocket  launch.  The  auroral 
clutter  continued  to  be  very  weak  and  patchy  since  their  reappearance 
around  0640  IT.  The  auroral  clutter  map  corresponding  to  approximately 
two  minutes  after  the  second  rocket  launch  is  shown  in  Figure  4.  The 
echoes  at  around  900  km  in  the  east  azimuths  may  be  remnants  of  the  0B 
seen  around  the  first  rocket  launch.  The  isolated  echo  patches  south 
of  Chatanika  and  seemingly  aligned  along  a magnetic  latitude  may  be 
related  to  a visual  arc,  although  all-sky-camera  photographs  available 
from  Chatanika  did  not  reveal  any  aurora  due  to  continued  cloud  coverage. 

2.  19  February  1971 — Nike  llvdac 

A Nike-Hydac  rocket  (Nil  74-1)  carrying  a TMA  chemical  release 

payload  and  a falling-sphere  atmospheric  density  experiment  was  launched 

at  0711  IT  during  quiet  conditions.  Geomagnetic  activity  at  College  was 

low,  although  some  perturbation  began  just  prior  to  launch  at  about 
25 

0720  FT.  The  magnetic  deflections  can  be  interpreted  in  terms  of  the 
devel>  jiment  of  a weak  westward  current,  located  poleward  of  College. 

This  current  system  was  localized,  since  Barrow  did  not  record  substantial 
geomagnetic  activity. 

The  auroral  clutter,  similar  to  those  seen  during  the  first 
rocket  launch,  formed  a weak  DU  approximately  200  km  northeast  of 
Chatanika  (and  directly  over  Fort  Vukon).  Shortly  before  the  rocket 
launch  (0742:30  UT),  the  I)U  bifurcated  into  two  latitudinal  segments, 
possibly  due  to  the  development  of  a visual  arc  within  the  DU.  (A  faint 
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auroral  arc  was  observed  on  the  northern  horizon  of  the  Chatanika  all- 

o 

skv-camera  photograph  and  also  by  a meridian-scanning  5577-A  photometer 

25 

operated  by  AFCKL/USF.  ) The  auroral  clutter  map  taken  approximately 
two  minutes  after  the  rocket  launch  is  shown  in  Figure  5.  By  this  time, 
the  echoes  had  weakened,  becoming  patchy. 

The  electron  densities  were  very  low,  as  measured  by  the 
25 

Chatanika  radar.  The  E-region  electron  density  was  approximately 
1 3 

10  el/cm  , again  leaving  open  the  possibility  that  the  electric  field 

could  have  been  greater  than  30  mV/m.  The  F-region  electron  density  was 

4 3 

no  greater  than  0.7  y 10  el/cm  , determined  only  alter  43  minutes  of 
integration . 

3.  21  February  197  1 — Black  Brant  (H1KIS) 

A Black  Brant  rocket  (A18. 116-1)  with  a HIR1S  payload  was 

2 

launched  at  0916:30  FT  during  active  auroral  and  geomagnetic  conditions. 
Unfortunately,  during  this  time  period  the  Homer  phased-array  radar 
developed  software  difficulties  associated  with  the  electronic  scanning. 
Auroral  clutter  data  were,  however,  collected  on  the  139-MHz  radar  but 
will  not  be  presented  due  to  difficulty  in  its  interpretation. 

1 . 25  February  1974— Black  Brant  (MULT I ) 

A Black  Brant  rocket  (A18. 219-1),  carrying  a MULTI  payload  to 

observe  upper-atmospheric  chemical  processes  and  IK  emissions  within  a 

bright  auroral  arc,  was  launched  at  0738:30  FT  during  disturbed  auroral 

25 

and  geomagnetic  conditions. 

The  characteristic  DB  radar  signature  similar  to  that  presented 
in  Figure  2 (Section  II)  was  seen  throughout  this  time  sector  around  the 
launch.  The  auroral  clutter  map  taken  approximately  two  minutes  after 
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the  rocket  launch  is  shown  in  Figure  6.  The  map  format  was  changed  from 
the  polar-plot  format  shown  in  the  first  three  figures  to  the  rectangular 
slant -versus-azimuth  format,  due  to  on-line  difficulties  with  the  polar- 
plot  format. 

The  1)1$  is  seen  to  extend  from  50  km  north  of  Chatanika  to 
approximately  200  km  south  of  Chatanika  where  the  auroral  clutter  is  cut 
oil  by  magnetic  aspect  sensitivity.  Examination  of  Chatanika  all-sky- 
camera  photographs  showed  that  a bright  visual  arc  was  located  along  the 
1)1$  poleward  boundary.  A second,  faint  arc  segment  was  also  seen  directly 
over  and  to  the  west  of  Chatanika  that  coincides  with  the  bifurcation  in 
the  I)B  over  Chatanika  (Figure  6).  The  bifurcation  was  apparently  produced 
by  the  depression  of  the  electric  field  by  particle  precipitation.  A 
westward-traveling  surge  passed  over  Chatanika  within  the  next  live 
minutes  along  the  I)L$  poleward  boundary. 

The  DB  is  typical  of  that  usually  observed  during  a moderately 
active  evening  period  when  the  1)1$  is  essentially  collocated  with  the 
eastward  electrojet  and  the  poleward  boundary  representing  the  lower 
latitudinal  boundary  of  the  Harang  discontinuity.  Furthermore,  the  DB 
represents  a region  of  downward  field-aligned  current  flow. 

In  accord  with  the  above  description,  the  Chatanika  radar 
measured  a northeastward-directed  electric  field  vector  with  a strength 

of  approximately  60  mV/m.  E-layer  electron  densities  were  of  the  order 

5 , 3 

of  10  el/cm  . 

1$.  Doppler  Capability 

A third  field  trip  was  made  during  November  1971  to  install  a hard- 
wired fast-Eourier-transform  (FFT)  processor  for  near  on-line  spectral 
observations  of  auroral  clutter.  The  computer  programming  necessary  to 
operate  the  cl  it ter-mapping  radar  and  to  take  data  in  a desired  format 
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FIGURE  6 REAL-TIME  B-SCAN  — 25  FEBRUARY  1974,  0740  46  UT 


was  completed  and  debugged.  The  programming  necessary  to  process  the 
prerecorded  clutter  data  through  the  KFT  processor  and  to  plot  Doppler- 
velocity  maps  was  also  completed.  In  the  process  ol  debugging  this 
latter  program,  it  was  found  that  the  FIT  processor  would  not  execute 
transforms  at  the  required  rate.  At  the  end  of  the  field  trip,  this 
processor  was  returned  to  the  SHI  Menlo  Park  facility.  Approximately 
12  hours  of  clutter  Doppler  data  were  recorded  during  the  field  trip. 

To  see  ll  the  data  had  been  recorded  properly,  a data  set  collected 
during  a period  of  substorm  activity  (1128  UT,  20  November  1971)  was 
processed  using  software  KFT  techniques.  The  results  are  shown  in 
Figure  7.  The  on-line  auroral  clutter  map  is  shown  in  Figure  7(a).  In 
order  to  accommodate  the  larger  computer  memory  requirements,  only  the 
data  from  400  km  to  1100  km  slant  range  are  processed  and  mapped.  The 
clutter  map  does  not  appear  to  reveal  any  complex  echo  structure  that 
might  be  expected  under  substorm  conditions.  Auroral  clutter  is  seen  to 
extend  from  450  km,  where  the  radar  echoes  are  cut  off  by  magnetic  aspect 
sensitivity,  to  1100  km,  where  the  clutter  appears  weaker  and  more  patchy. 

o 

To  minimize  computation  time,  only  the  data  taken  every  2 m eleva- 
o 

t ion  and  4 in  azimuth  were  processed.  Contours  of  constant  Doppler 
velocity  were  then  hand-drawn  through  the  mean  Doppler-veiocity  values. 

The  results,  shown  in  Figure  7(b),  are  clearly  more  complex  in  velocity 
space  than  they  are  in  clutter-amplitude  space.  There  appears  to  be  no 
obvious  relationship  between  the  observed  radial  Doppler  velocities  and 
clutter  amplitude. 

Other  than  the  fact  that  the  data  were  indeed  recorded  properly, 
little  physics  can  be  extracted  from  this  particular  example.  Actual 
analysis  and  utility  of  the  Doppler  data  will  have  to  begin  with  more 
stable  periods — e.g.,  those  in  which  there  exists  a well-behaved  after- 
noon or  evening  DB  or  stable  discrete  arcs  (DA). 
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AURORAL  CLUTTER  MAP  AND  CORRESPONDING  DOPPLER 
VELOCITY  MAP  — 1128  UT,  20  NOVEMBER  1974 


St-arch  !<>r  Plasma  Pine  Backseat  ter 


C. 

A linal  1 ielcl  trip  was  made  during  10-17  March  1975  to  search  for 

the  possible  existence  of  plasma-line  backscatter  during  periods  of 

intense  auroral  activity.  This  study  was  initiated  by  Dr.  John  Minkoff 

ol  Hiverside  Research  Institute.  Having  participated  in  ionospheric  RF 

27 

heating  experiments  where  plasma-line  backscatter  was  observed,  Dr. 
Minkoff  approached  the  National  Science  Foundation  with  a proposal,  which 
was  accepted,  to  search  for  the  possible  presence  of  plasma-line  back- 
scatter that  might  be  associated  with  auroral  activity. 

The  139-  and  398-MHz  radars  at  Homer  were  operated  using  the  60-ft 
reflector  antenna.  A Hewlett-Packard  spectrum  analyzer  was  connected 
to  the  30-MHz  IF  of  one  of  the  two  radars.  The  analyzer  allowed  a search 
in  frequency  on  both  sides  of  the  "on  frequency”  auroral  clutter  echo. 
Searches  were  made  to  +10  MHz,  A range  gate  was  used  to  allow  only  the 
auroral  clutter  into  the  spectrum  analyzer. 

The  operational  procedure  was  to  make  a search  for  auroral  clutter 
by  manually  scanning  the  antenna  and  using  the  normal  A-scan  display. 
After  a large  signal-to-noise  ratio  (SNR)  echo  was  found,  the  manually 
controlled  range  gate  was  positioned  on  the  auroral  echo  and  the  spectrum 
analyzer  was  used  at  both  139  and  398  MHz,  alternately. 

Over  a period  of  six  days  this  technique  was  used  to  search  for 
plasma-line  backscatter  at  all  hours  of  the  night  and  day.  During  the 
six  days,  both  diffuse  and  discrete  auroral  echoes  were  observed  with 
amplitudes  up  to  40  dB  SNR.  These  signal  strengths  correspond  to  some 
of  the  strongest  auroral  clutter  ever  observed  with  the  Homer  radar. 
However,  at  no  time  were  echoes  observed  off-channel  at  a frequency  that 
might  correspond  to  the  plasma  frequency.  The  search  included  not  only 
the  E layer  with  auroral  clutter  present,  but  the  F layer  above  existing 
auroral  clutter. 


IflBIW-Y* 


The  negative  results  suggest  at  least  two  possibilities.  First, 
if  we  assume  that  plasma-line  backscatter  is  related  to  the  amount  of 
energy  that  is  deposited,  it  can  be  shown  through  thatanika  radar  measure- 
ments that  Joule  and  particle  heating  are  easily  several  orders  of  magni- 
tude less  than  that  produced  by  the  RF  heater  used  in  the  coherent  heating 
experiments.  And  second,  the  tact  that  the  aurora  constitutes  an 
incoherent  source  (i.e.,  broadband  in  frequency  space)  of  energy  suggests 
that  sideband  generation  at  a particular  frequency — say,  the  plasma 
frequency — does  not  take  place.  Instead,  any  sideband  generation,  if 
present,  may  have  a broadband  distribution,  similar  to  the  input  energy 
distribution,  that  is  not  intense  enough  to  be  detected  at  any  particular 
frequency  by  an  auroral  backscatter  radar.  In  this  case,  integration  in 
frequency  space  (as  well  as  in  time)  excluding  the  "on  frequency"  returns 
may  reveal  a detectable  signal. 
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IV  Till  ION-CYCLOTRON  WAVE  AND  ASSOCIATED  INSTABILITIES 


A.  General 
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Recently  D'Angelo  suggested  that  the  Type  III  spectra  associated 
29 

with  the  radar  aurora  may  possibly  be  explained  by  unstable  ion  cyclo- 
tron waves.  A number  of  previous  works  on  the  electron  cyclotron  waves 
indicate  that  field-aligned  streams  and  currents  can  make  the  wave  grow. 
It  is  expected  that  similar  instabilities  are  also  possible  for  the  ion 
cyclotron  waves.  Recently,  a large  amount  of  experimental  evidence  has 
been  presented  confirming  the  presence  of  field-aligned  current  and 
streams  in  the  auroral  region.  It  is  therefore  important  to  investigate 
possible  instabilities  associated  with  ion  cyclotron  waves. 


We  have  started  to  work  on  this  problem  during  the  current  contract, 
but  the  work  is  not  yet  completed.  At  the  present  time,  various  forms 
of  dispersion  equations  based  on  kinetic  theory  have  been  derived,  and 
the  propagation  cha racteris t ics  for  the  ion  cyclotron  waves  in  a stable 
plasma  have  been  determined.  These  results  are  presented  in  this  section. 

B.  Dispersion  Equation 

In  a uniform,  collisionless  magnetoplasma,  the  dispersion  equation 

— • — « 

for  a quasi-static  wave  (E||k)  can  be  derived  from  the  Vlasov  and  Poisson 
equations.  The  procedures  are  as  follows:  One  first  integrates  the 

Vlasov  equation  along  the  orbit  of  the  particles  to  obtain  a charge 
expression  for  each  species  in  terms  of  the  wave  field.  Then  the  charge 
expression  is  substituted  into  Poisson's  equation  to  eliminate  the  wave 
field  and  to  obtain  a dispersion  equation.  By  assuming  that  only  first- 
order  perturbations  are  important,  the  Vlasov  equation  and  Poisson's 
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equation,  as  well  as  the  above  derivation  procedure,  can  be  linearized. 

By  this  we  mean  that  the  perturbing  wave  field  K is  a first-order  (small) 

quantity,  and  the  particle  distribution  function  can  be  written  as  f - 

f + f , where  ( is  the  equilibrium  distribution,  and  f is  the  first- 
0 1 0 _ 1 

order  perturbed  distribution  due  to  E.  In  the  integration  of  the 
linearized  Vlasov  equation,  the  integration  is  performed  along  the  orbit 
of  the  unperturbed  particles. 


While  the  procedure  is  as  described  above,  we  omit  the  lengthy 


general  form  of  the  dispersion  equation: 
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where 

^ = Plasma  frequency 

P 

w Gy ro frequency 
c 

k = Wave  number 

J (*)  = n -order  Bessel  function  of  first  kind 
n 

f - Unperturbed  particle  distribution  function 
X = kxVxA> 

c 

n = Integer. 

The  subscripts  " i,  and  "j."  stand  for  the  components  parallel  to  and 

perpendicular  to  the  static  magnetic  field  B , respectively.  £ stands 

o s 

for  summing  over  all  charged-particle  species. 

For  the  case  of  collisions,  the  Boltzman  equation  must  be  used  to 
replace  the  Vlasov  equation.  However,  the  exact  integral  form  of  the  rate 
of  change  of  the  distribution  function  due  to  the  collisions  is  extremely 
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complicated.  A usual  approximation  to  that  integral  form  is  known  as  the 
IK  I K model,  where  the  collisional  integral  in  the  Boltzmann  equation  is 
replaced  by  - f()(l  + n^/n^)  - f].  Here,  n^  is  the  unperturbed  particle 
density,  n ^ is  the  first-order  perturbed  density,  and  ' is  the  effective 
collision  frequency.  The  BGK  collision  model  allows  the  particles  to 
relax  in  position  space  to  the  local  density,  rather  than  to  the  unper- 
turbed density  n^  which  corresponds  to  f ^ , and  conserves  particles  in 
position  space. 

With  the  BGK  model  of  collisions  lor  each  species  taken  into  account, 

the  derivation  of  the  dispersion  equation  is  similar  to  that  for  the 

collisionless  case.  It  can  be  shown  that  the  general  form  of  the 

2 

dispersion  equation  is  again  given  by  Kq.  (1)  after  replacing  w by 

[l  - j(V/<u)] 

P 

and  replacing  a in  - k ' V , - n*  ) by  (<*>  - jv),  where  j -1  . 

c 

While  the  general  form  of  the  dispersion  equation  is  given  above, 
we  now  proceed  to  reduce  Kq.  (1)  for  a number  of  specific  unperturbed 
distribution  functions. 

1 . Zero-Drift,  Isotropic  Maxwellian  Distributions 

The  distribution  functions  for  the  electrons  and  ions  take 

the  form  ,2  2v 

(Yx+vi; ) 


where  the  subscripts  for  the  electrons  and  ions  are  omitted  and  where 
stands  for  the  thermal  speed. 
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Substituting  F.q.  (2)  into  Eq.  (1)  and  performing  the  integrations 


in  velocity  space  yields  the  dispersion  equation: 

2 -X 


£ kl  — I 

e , i 2 n--° 

(JU 

c 


1 (X) 

n 


1 


^ k»Ut 


ZC  ) 


= 0 


(3) 


2 ..  th 

where  ' ( k if  /x  ) , I (X)  n -order  Bessel  function  ol  the  third 

t c n 
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kind,  and  Z(i)  is  known  as  the  plasma  dispersion  lunction, 

! - (x-n  x )/,2  k,|U  . 

c t 

When  the  collisions  are  included,  Eq.  (3)  becomes 


2 v -X  , 

x (1-j  -)  e I (X)  f 

2 r 2 p x S n 1 

k + £ ki.  — £ : 

e , i 2 n=-°°  X 


1 + — 


^ k"Ut 


ZC) 


= 0 


(4) 


where  1 in  Eq.  (!)  now  becomes 

u)  - jv  - n1*) 


■J 2 k ||U 


2.  Anisotropic,  Zero-Urift  Maxwellian  Distributions 


In  this  case  the  thermal  velocity  takes  different  component 
values  in  the  directions  parallel  to  and  perpendicular  to  the  static 
magnetic  field.  The  Maxwellian  distributions  for  the  electrons  and  ions 
take  the  following  form:  2 

vj.  vji 


f (V) 
o 


2 2 
1 2Ui  2U  || 

3/2  2 e 

(2n)  U||Ux 


(5) 


where  l1 1| , i stands  for  the  thermal  velocity  parallel  and  perpendicular  to 

B , respectively, 
o 

Inserting  Eq.  (5)  into  Eq.  (1)  yields 
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2 ‘> 
k~  + 2 ki 

e,  i 


2 2 -X 

m l'i  cl  (X) 

P v-  n 


2 2 n=-“ 

u(1  -c 


U || 

w - (1 £)  n«> 

2 c 
Uj.  c 

1 -t  Z(s) 

-.2  k ii  U || 
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where  X and  j are  now  defined  as 


k±lri 

ID 
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oj  - nui 


v 2 k U II 


In  the  limit  of  Ux  = U||  = , Eq.  (5)  reduces  to  Eq.  (2).  It 

is  clear  that  Eq.  (6)  becomes  identical  to  Eq.  (3). 

2 

For  the  case  of  collisions  included,  we  simply  replace  by 

2 v P 

u>  (l  - )— ) and  '•*>  by  (a,-jv),  in  Eq.  (6)  to  obtain  the  new  dispersion 

p 

equation.  In  comparison  with  the  dispersion  equation  derived  by  Tataronis 
and  Crawford,  1 we  point  out  that  Eq.  (19)  of  their  paper  appears  to  be 
incorrect . 


3.  Field-Aligned  Drift  Maxwell  tans 


Let  U be  the  drift  velocity  along  B ; the  distribution  function 
d o 


is  given  by 


[Vx  + (V  ||  - ud>2] 


Vv) 


3/2 


2U 


(7) 


and  the  dispersion  equation  lEq.  (3)]  becomes 
2 

e I (A)  | 

® n I 
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uj  e I (X) 
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where  i*>  - <*>  - k lT 


the  Doppler-shifted  frequency,  and 


UJ  — fjUU 


V2  k||U 


Similarly,  if  V||  in  Eq.  (5)  is  replaced  by  V ||  - U , the  corresponding 

d 

dispersion  equation  is  again  given  by  Eq.  (6),  with  replaced  by 

u/  - u)  - k||U  . 


1 . Ambient  Plasma  Plus  a Field-Aligned  Stream 

Let  D (lD,k)  - 0 be  the  dispersion  equation  for  the  ambient 
P 

plasma — i.e.,  let  D (u^k)  represent  the  left-hand  side  of  any  dispersion 

— * 

equation  in  the  previous  sections.  If  we  let  D (u>,k)  be  the  dispersion 

s 

function  due  to  the  stream,  the  overall  dispersion  equation  of  the 
system  is  then  given  by 


D (w,k)  + D (w,k)  = 0 ( 

P s 

To  derive  the  expression  D ( '-u , k ) , we  assume  that  the  distri- 

s 

bution  of  the  stream  takes  the  form: 


f (V)  = fj(Vi.)  6 (V | - U 
s s 


where  U is  the  stream  velocity  along  B . Substituting  Eq.  (10)  into 
s o 

the  integral  of  Eq.  (1)  yields 

nil  2 

v - r-  s k|, 

D (w,k)  = u>  £ I 11 I 

s s n=-®  u>-k||U  -n \i  1 u)-k||U  -n  2 

s s s s 


where  a)  = stream  plasma  frequency,  (carrying  the  sign  of  charge)  = 
s s 

gy rof requency  of  the  particles  in  the  stream,  and  I ^ and  I are  given  by 


a 2 dfj.  2 

I = 2ix  I J (X  ) dVj  ; 1=2*  J^(X  ) fj.Vj.dVj. 

1 _ n s dVx  2 ".ns 


For  the  case  ol  a munoenerKtt ic  stream  along  the  field  line,  f takes  the 

s 

form 


f 

s 


1 

2 X l 


5(VX)  (V  -U  ) 
b 


(13) 


Integration  of  1^,  I ^ in  Eq.  (12)  is  then  straightforward.  Equation  (11) 
can  be  reduced  to  the  form: 


I)  C,k) 

s 


(1  1) 


C . Ion  Cyclotron  IVaves 

In  this  section  we  examine  the  propagation  characteristics  of  ion 
cyclotron  waves  in  a stable  plasma.  For  this  purpose  we  assume  that  the 
electrons  and  ions  exhibit  an  isotropic,  stationary  Maxwellian  distri- 
bution, and  we  use  the  dispersion  equation  given  by  Eq.  (3)  or  Eq.  (4). 

1 . t o I 1 i s ion  1 ess  Case 

It)  t ho  case  ol  no  collisions,  the  dispersion  equation  is  given 
by  1 q . (3).  For  the  isotropic,  stationary  Maxwellian  plasma,  the  wave 
o.,«  - < mtained  in  t hi'  dispersion  equation  consist  of  no  instabilities. 

I hi  *avi  may  be  subjected  to  I-andau  damping.  In  the  linear  regime,  I.andau 
damping  is  proportional  to  ImZ(l).  (Im  stands  for  the  imaginary  part.) 

For  thi  wave  propagating  at  a large  angle  with  respect  to  B (k|i  becomes 
small),  lm/("  ) is  then  much  smaller  than  KcZ(L),  and  the  waves  are  sub- 
stantially free  of  I^indau  damping.  In  the  remainder  of  this  section  we 
shall  discuss  the  three  cases:  k||  = 0;  kj.  0;  and  kj  ^ 0,  k||  t 0,  and 

kj.  » kn. 

32 


0 


For  this  case,  the  propagation  direction  is  exactly  per- 

oendicular  to  U (kiB  ),  In  the  following,  we  use  the  subscripts  "e" 
o o 

and  "i"  to  designate  the  parameters  for  the  electrons  and  ions,  respect- 
ively. Solving  for  ion  cyclotron  waves,  we  assume  that 


m T 

e e _e 

m T 
i i i 


which  yields  the  approximation 


1 < a ) exp  - X ~ (1  - X ) 

n e n e n 

2 n ! 


Using  F.qs . (15)  and  (16)  we  can  approximate  Eq.  (3)  by  the  form: 

2 22  T . 
u:  2n  u)  I ( X ) 

i - |,e-  - - i Si H—L-  exp  . x = 0 . (] 

2 2 n=l  2 2 X l 

tx  - tx  id  - (nui  ) t 

ce  ci 

t h 2 2 

For  the  n ion  cyclotron  harmonic  wave,  yields 

2 2 2 

2n  ^ I (X  ) 

2 2 pi  ce  n l . , , 

<x  ( n<t>  ) + — r ex p - X _ U 

ci  2 A t 

x i 
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2 2 2 

where  tx  = <x  + ta  and  ntx  < uj  < (n+l)u> 

u pe  ce  c i ci 


Equation  (18)  can  be  plotted  to  yield  the  dispersion  curves  — i.e.,  plots 
of  (x  versus  kj.  for  each  n and  the  other  plasma  parameters.  From  the  x_ 
versus-kj.  plots,  the  phase  velocity  tx/k^  can  be  determined.  Two  limiting 
approximations  can  be  made  to  simplify  E;q.  (19). 


i 


Approximation  1:  « 1.  In  the  case  of  phase  velocity 

much  larger  than  the  ion  thermal  speed  U^,  or  wavelength  much  larger  than 
the  ion  I.armor  radius,  we  have  ^ « 1.  Equation  (18)  can  be  approximated 

l 

by 
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2 2 2 
n u)  id 

2 .2  pi  ce 

ut  (n’x  ) + 

c i 2 
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Approximation  2:  ^ » 1.  Reversing  the  limiting  case 

discussed  above,  we  use  1 (x)  ' (exp  x)/.  2tx.  Equation  (18)  can  then  be 

n 

approximated  by 
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(20) 


b.  kx  = 0 


magne  t ic 


For  the  case  of 
field,  the  dispersion 


wave  propagation  exactly  along  the  static 
equation,  Eq.  (3),  reduces  to 


ID 
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(21) 


where  t'  (!)  = dz(;,)/d^. 

It  is  seen  that  the  cyclotron  frequencies  have  disappeared 
from  the  dispersion  equation.  The  solutions  of  Eq.  (21)  are  simply  the 
longitudinal  plasma  waves,  which  degenerate  into  plasma  oscillation  when 
Te  Ti  = 0.  Since  the  imaginary  part  of  7.'  is  quite  large,  these  plasma 
waves  are  subject  to  strong  Ixtndau  damping. 


J 
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kj.  / 0,  k , / 0,  and  kj  k 


c . 

In  this  case,  we  consider  that  ion  cyclotron  waves  propa- 
gate with  la:  ye  angle  witli  respect  to  B . Since  kr  is  small,  1 and  ' 

o e i 

are  much  greater  than  unity.  Thus  the  Eandau  damping  can  again  be 
neglected.  After  some  algebra,  the  dispersion  equation  can  be  reduced  to 
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(22) 


where  0 is  polar  angle  between  B and  k. 

o 

th 


For  n ion  cyclotron  harmonic  waves,  Eq.  (22)  yields 
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where  n<D  < id  < (n+  1 ) id 

ci  ci 


The  asympotic  forms  for  A « 1 and  A » 1 are  similar  to  that  for  the 

i i 

case  of  kj  = 0,  with  the  factor  sin  0 multiplying  id  and  id  as 

ce  pe 

indicated  in  Eq.  ('22). 


2 . Finite  Collisions 

When  the  effects  of  collisions  cannot  be  neglected,  the  ion 

cyclotron  waves  are  subject  to  collisional  damping.  The  dispersion 

equation  for  ion  cyclotron  modes  propagating  perpendicular  to  B is  now 

o 

given  by 
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Future  Work 


\ 


In  previous  sections  we  have  presented  a number  of  general  forms 
of  the  dispersion  equation  under  the  various  plasma  conditions.  For  the 
case  of  a stable  plasma,  the  dispersion  equation  has  been  solved  for  the 
ion-cvclotron  harmonic  waves.  When  the  finite  collisions  are  included, 
the  formulas  of  the  collisional  damping  rate  for  the  ion  cyclotron 
harmonic  waves  have  also  been  derived. 

As  the  next  step,  we  should  consider  several  models  of  unstable 
plasma  in  the  auroral  region,  and  investigate  the  instabilities  of  the 
ion  cyclotron  waves.  With  and  without  collisions,  the  onset  (threshold) 
conditions  for  instability  and  the  formulas  of  the  growth  rate  should 
be  derived. 

It  is  expected  that  the  above-mentioned  quantities  of  the  ion 
cyclotron  waves  are  closely  related  to  the  field-aligned  streams  and 
currents.  Thus,  these  results  can  be  used  not  only  to  explain  the  Doppler 
spectra  of  auroral  clutter,  but  also  as  a diagnostic  tool  to  measure 
the  field-aligned  currents. 


1 


37 


f 


I 


V DISCUSSION 


It  is  clear  from  the  results  presented  in  this  report  that  great 
strides  have  been  made  during  the  past  2 or  3 years  toward  an  under- 
standing of  auroral-clutter  production  mechanisms.  With  the  instal- 
lation of  both  the  398-MHz  phased-array  radar  at  Homer  and  the  incoherent- 
scatter  radar  at  Chatanika,  it  became  possible  for  the  first  time  to 
continuously  monitor  the  plasma  parameters  within  a radar  scattering 
volume  producing  auroral  clutter.  The  combined  use  of  the  two  radars 
together  with  other  auroral/magnetospheric  measurements  has  made  it 
possible  to  evaluate  existing  plasma  instability  models  as  candidate 
scattering  mechanisms  and  to  determine  auroral  clutter  characteristics 
that  must  be  explained  by  any  successful  theory. 

Recent  results  showed  that  the  Buneman-Farley  two-stream  instability 

is  still  a prime  candidate  as  the  source  of  auroral  clutter.  However, 

it  is  now  clear  that  while  the  two-stream  primary  waves  may  in  fact  be 

the  driving  force,  auroral  clutter  is  apparently  not  produced  by  direct 

backseat  ter  from  these  irregularities.  Evidence  for  the  two-stream 

instability  as  the  prime  driving  source  is:  (1)  the  existence  of  a 

15 

threshold  electric  field  strength  oi  30  mV/m;  and  (2)  Doppler  velocity 

7 

characteristics  us  related  to  the  ion-acoustic  speed,  the  magnetic 

7,32  6,7 

aspect  angle,  and  the  current  direction. 

The  gradient-drift  instability  also  remains  a candidate,  however, 
primarily  through  its  role  in  the  production  of  secondary  plasma  waves. 

In  its  linear  form  this  instability  appears  incapable  of  producing 
scale  sizes  necessary  to  explain  auroral  clutter  observed  at  frequencies 
as  high  as  3 GHz.  The  only  model  that  allows  for  the  generation  of  small- 
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size  grad ient-dr i f t waves  is  still  not  completely  satisfactory.  As 

6 

discussed  by  Tsunoda,  the  ion  velocity  effects  that  appear  in  observed 

Doppler- veloc i ty  data  cannot  be  accounted  for  by  this  theory.  Furthermore, 

the  necessary  assumption  that  wave  steepening  does  not  occur  appears  to 

7 

break  down  at  large  primary  wavenumbers. 

The  question  of  whether  the  ion-cyclotron  instability  might  play 
a role  in  the  production  of  auroral  clutter  or  at  least  contribute  to 
the  modification  of  the  observed  echo  characteristics  remains  to  be 
investigated.  The  discovery  that  the  diffuse  auroral  clutter  (DB)  found 
in  the  evening  sector  is  collocated  with  downward  field-aligned  currents 
suggests  the  possibility  that  ion-cyclotron  waves  may  be  present. 

However,  the  fact  that  most  of  the  downward  field-aligned  currents  in 
the  DB  region  are  probably  carried  by  upward-moving  thermal  electrons  in 
the  F region  may  negate  this  suggestion.  There  does  remain  the  possibility 
that  visual  aurora  and  discrete  auroral  clutter  are  associated  with  ion- 
cyclotron  waves. 

The  determination  of  microscale  and  macroscale  relationships  (see 
Section  II)  has  also  converted  the  backseat  ter  radar  into  a valuable 
diagnostic  tool  for  auroral  and  magnetospher ic  research. 

The  typical  morphological  features  that  can  be  extracted  from  the 
auroral  clutter  were  presented  in  Section  II  (see  Figure  2).  Regions  of 
enhanced  electric  fields  and  significant  electron  density  are  mapped 
by  the  auroral  clutter.  Both  Doppler  and  amplitude  characteristics  of 
auroral  clutter  appear  capable  of  providing  quantitative  information 
regarding  the  electric  field  and  E-region  current  density. 

A direct  application  of  these  morphological  features  is  to  test 
the  NRl./DNA  Auroral  codes.  The  capabilities  of  the  Homer  (609)  and 
Chatanika  (617)  radars  have  been  extended  so  that  they  are  able  to  monitor 
nearly  all  of  the  ionospheric  plasma  parameters  that  ore  the  input  and 
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output  ol  the  NK1./DNA  Auroral  codes.  The  only  key  parameter  *ot  presently 
measurable  by  the  ground-based  radars  is  the  1 ield-al igned  current,  which 
is  monitored  by  a vector  magnetometer  on  the  Navy  TRIAD  satellite.  (At 
least  two  other  satellites  are  able  to  provide  similar  measurements — 
ISIS-2  and  the  Air  Force  S3-2.) 

The  principal  outputs  of  the  Auroral  code  are  the  magnitude  and 
spatial  variation  ol  (1)  the  ionospheric  electric  field,  (2)  the  electro- 
jet  currents,  and  (3)  the  net  electric  current  flow  parallel  to  the 
geomagnetic  field  line  (field-aligned  current).  The  617  and  609  radars 
can  measure  the  ionospheric  electric  field  and  electrojet  currents  as  a 
function  of  latitude.  At  the  present  time  it  is  not  possible  to  measure 
field-aligned  current  with  the  ground-based  systems. 

The  principal  role  of  the  Homer  radar  in  such  a coordinated  experi- 
ment is  illustrated  in  Figure  8.  The  Chatanika  radar  measurements  are 
generally  restricted  to  the  region  in  or  near  the  local  geomagnetic 
meridian.  The  plane  in  which  data  are  gathered  by  a satellite  rarely, 
if  ever,  is  coincident  with  the  Chatanika  meridian.  Consequently,  the 
two-dimensional  (latitude,  longitude)  mapping  capability  of  the  609  radar 
is  the  essential  link  between  the  617  radar  and  TRIAD  measurements.  The 
sensitivity  of  the  Homer  radar  to  regions  of  enhanced  electric  field 
strengths  and  electrojet  currents  can  be  used  to  map  these  regions, 
allowing  accurate  spatial  (and  temporal)  extrapolat ion  between  the  lanes 
of  the  TRIAD  pass  and  the  Chatanika  radar  measurements.  Note  in  Figure  8 

that  by  utilization  of  the  mechanical  steering  capability  of  the  Homer 

o 

radar  in  addition  to  its  10  electronic  scan,  the  longitudinal  coverage 
can  be  significantly  extended. 

In  the  course  of  the  reported  auroral  clutter  studies,  related 

studies  have  shown  that  auroral  clutter  is  probably  not  limited  only  to 

33 

the  auroral  E region.  Ott  and  Farley  recently  developed  a plasma 
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FIGURE  8 A TECHNIQUE  FOR  USING  THE  HOMER  PHASED-ARRAY  RADAR  TO  ALLOW  ACCURATE  EXTRAPOLATION 
BETWEEN  SATELLITE  MEASUREMENTS  AND  CHATANIKA  RADAR  MEASUREMENTS 


instability  model  that  predicts  the  generation  of  field-aligned  irregu- 
larities in  the  auroral  F region.  The  instability  is  based  on  an  aniso- 
tropic ion  velocity  distribution  brought  about  by  charge-exchange  processes 
in  the  F region.  The  driving  force  is  again  a large  electric  field 
(’15  mV/m).  The  discovery  of  this  instability  is  a prime  example 
illustrating  our  still  rather  meager  understanding  of  electric  field 
effects  and  plasma  instabilities. 

The  irregularities  associated  with  this  instability  have  apparently 

gone  undetected  by  most  auroral  backscatter  radars  due  to  its  field- 

aligned  characteristics.  Because  of  the  geometry  of  the  geomagnetic 

field  at  high  latitudes,  it  is  virtually  impossible  for  ground-based 

radars  to  view  the  F region  in  directions  nearly  perpendicular  to  the 

field  lines.  However,  there  is  experimental  evidence  that  seems  to 
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confirm  its  existence.  I<eadabrand  et  al.  made  auroral  radar  measure- 
ments at  frequencies  of  lO]  and  800  MHz  that  showed  echoes  up  to  an 
altitude  of  160  km  or  so,  which  was  the  maximum  height  at  which 

orthogonality  was  achieved.  A second,  non-auroral  radar  observation 
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was  made  at  106.1  MHz  from  Stanford,  California.  During  periods  of 
high  magnetic  activity  when  the  auroral  oval  shifted  to  low  latitudes, 
field-aligned  echoes  were  observed  to  altitudes  as  high  as  300  km. 

Although  the  instability  proposed  by  Ott  and  Farley  does  not  predict 
backscatter  as  low  as  106  MHz,  the  observation  clearly  demonstrates 
that  F-region  auroral  clutter  does  indeed  occur. 
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VI  CONC IUS IONS  ANI)  RECOMMENDATIONS 

We  have  shown  that  auroral-clutter-producing  mechanisms  can  be 
effectively  studied  through  coordinated  field  experiments  utilizing  the 
Homer  609  radar,  the  Chatanika  617  radar,  and  other  supporting  measure- 
ments. Most  of  the  results  obtained  to  date  pertained  to  diffuse  auroral 
clutter.  The  choice  of  the  diffuse  auroral  clutter  for  the  initial 
study  was  based  on  its  well-behaved  spatial  and  temporal  characteristics. 
In  contrast,  the  discrete  auroral  clutter  is  characterized  by  a more 
complex  and  dynamic  behavior,  often  associated  with  auroral  substorm 
activity. 

We  have  shown  that  difluse  auroral  clutter  is  characterized  by  a 
threshold  electric  field  strength  of  30  mV/m.  Diffuse  auroral  clutter 
was  also  shown  to  be  collocated  with  regions  of  diffuse  particle  preci- 
pitation and  of  auroral  electrojet  flow.  Evidence  was  also  presented  to 
indicate  that  the  amplitude  is  probably  proportional  to  the  current 

density,  in  a manner  similar  to  that  observed  at  50  Mllz  by  Greenwald  et 
, 17,18 

al.  Both  amplitude  and  Doppler  characteristics  were  found  to  be 

generally  consistent  with  the  production  of  field-aligned  irregularities 

13 

bv  a plasma  instability  model  not  unlike  that  presented  by  Sudan  et  al. 

The  discrete  auroral  clutter  remains  to  be  characterized.  With  the 
phased-array  radar,  it  will  be  possible  for  the  first  time  to  characterize 
the  discrete  auroral  clutter,  in  both  amplitude  and  Doppler.  The 
clutter-producing  mechanism  associated  with  the  discrete  echoes  are 
believed  to  be  similar  to  the  field-aligned  nuclear  clutter  observed  due 
to  beta  ionization. 
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The  new  theoretical  evidence  that  field-aligned  F-region  clutter 
can  be  expected  when  electric  field  strengths  of  the  order  of  50  mV/m 
or  greater  exist  opens  a new  area  of  study  pertinent  to  clutter  environ- 
ments. The  theoretical  results  appear  to  be  supported  by  at  least  two 
past  experimental  studies  conducted  by  Stanford  Research  Institute. 

Finally,  advantage  should  be  taken  of  the  Homer  and  Chatanika 
radars  to  test  the  NRL/DNA  Auroral  codes.  A discussion  on  how  the 
Homer  and  Chatanika  radars  could  be  utilized  to  test  the  NRL./DNA 
Auroral  codes  was  presented  in  Section  V. 
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ATTN:  Conrad  L.  Longmire 
ATTN:  Dave  Sowfe 
ATTN:  Ralph  KUb 
ATTN:  Gary  McCarter 
ATTN:  P.  Fischer 
ATTN:  M.  Schcibe 

The  Mitre  Corporation 
AT'  Bill  Mehuron 

North  Carolina  State  University  at  Raleigh 
ATTN:  Sec.  Officer  for  Walter  A . Flood 

Photometries,  Inc. 

ATTN:  Irving  L.  Kofsky 

Physical  Dynamics,  Inc. 

ATTN:  Joseph  B.  Workman 
ATTN:  A.  Thompson 

Physical  Sciences,  Inc. 

ATTN:  Kurt  Wray 

The  Trustees  of  Princeton  University 

ATTN:  F.  W Perkins,  Plasma  Physics  Lab 

R & D Associates 

ATTN:  R.P.Turco 
ATTN:  Forest  Gilmore 
ATTN:  H.  A.  Ory 
ATTN;  Robe rt  E.  LeLevier 
ATTN:  Richard  Latter 

R & D Associates 

ATTN:  Herbert  J.  Mitchell 

The  Rand  Corporation 

ATTN:  James  Oakley 
ATTN;  Cullen  Crain 

Raytheon  Company 

ATTN:  G.  D.  Thome 

RCA  Corporation 

Government  & Commercial  Systems 
ATTN:  lech.  Lib. 

ATTN  R.  D.  Bachinshy 

Riverside  Research  Institute 
ATTN:  R.  Popolow 

Science  Applications,  Inc. 

ATTN:  Robert  W.  Lowen 
ATTN:  Daniel  A.  Hamlin 
ATTN;  Lewis  M Linson 
ATTN:  D.  Sachs 

Space  Data  Co rpo ration 

ATTN:  Edward  F.  Allen 

Stanford  Research  Institute 

ATTN:  Warren  W.  Beming 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


Stanford  Research  Institute 

ATTN:  Ray  L.  Leadabrand 
ATTN:  M.  Baron 
ATTN:  Walter  G.  Chcsnut 
ATTN:  Robert  S.  Leonard 
ATTN:  F.  J.  Fremouw 
ATTN:  L.  L.  Cobb 
ATTN:  R.T.  Tsunoda 
ATTN;  R.  !.  P re  snell 
ATTN:  T.  N.  C.  Wang 
50  c>  ATTN:  Personnel 

Stanford  Research  Institute 
ATTN:  Dale  H.  Divis 

Technology  International  Corporation 
ATTN:  W.  P.  Hoquist 


Thiokol  Chemical  Corp.  Astro  Met  Site 
ATTN:  G.  C.  Alford 

Utah  State  University 

ATTN:  IX) ran  Baker 
ATTN:  C.  Wyatt 
ATTN : Kay  Baker 
ATTN:  D.  Burt 

Visidyne,  Inc. 

ATTN:  William  Reidy 
ATTN:  T.  C.  Degges 
ATTN : ,!  W . C arpentr  r 
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